Gluon Shadowing and Hadron Production at RHIC 
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Hadron multiplicity in the central rapidity region of high-energy heavy-ion collisions is investi- 
gated within a two- component mini-jet model which consists of soft and semi-hard particle produc- 
tion. The hard contribution from mini-jets is reevaluated using the latest parameterization of parton 
distributions and nuclear shadowing. The energy dependence of the experimental data from RHIC 
requires a strong nuclear shadowing of the gluon distribution in this model. The centrality depen- 
dence of the hadron multiplicity at ^/s = 130 GeV is reproduced well with the impact-parameter 
dependent parton shadowing. However, energy variation of the centrality dependence is needed to 
distinguish different particle production mechanisms such as the parton saturation model. 
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Formation of quark-gluon plasma (QGP) in high- 
energy heavy-ion coUisions hinges crucially on the initial 
condition that is reached in the earliest stage of the vio- 
lent nuclear interaction. Though many proposed signals 
can provide more direct measurements of the initial par- 
ton density, they must compliment results inferred indi- 
rectly from the measurement of final hadron multiplicity 
using either simple scenarios such as the Bjorken model 
101 or other dynamic models. Therefore, global observ- 
ables such as the rapidity density of hadron multiplicity 
can provide an important link of a puzzle that can eventu- 
ally lead one to a more complete picture of the dynamics 
of heavy-ion collisions and formation of QGP. Further- 
more, the study of energyand centrality dependence of 
central rapidity density (2) can also provide important 
constraints on models of initial entropy production and 
shed lights on the initial parton distributions in nuclei. 
For example, the available RHIC experimental data 
can already rule out the simple two-component model 
without nuclear modification of the parton distributions 
in nuclei In this paper, we will study within the 
two-component model how the RHIC data constrain the 
unknown nuclear shadowing of the gluon distribution in 
nuclei and how to further distinguish such a conventional 
parton production mechanism from other novel physics 
such as parton saturation 

Mini-jet production in a two-component model has 
long been proposed to explain the energy dependence 
of total cross section and particle production 

lOI,^ in high -energ y hadron collisions. It has also 
and incorporated in the HIJING 
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been proposed |15 

model 0,11^ to describe initial parton production in 
high-energy heavy-ion collisions. In this simple two- 
component model, one assumes that events of nucleon- 
nucleon collisions at high energy can be divided into those 
with and without hard or semi-hard processes of jet pro- 
duction. The soft and hard processes are separated by a 
cut-off scale pq. While the cross section of soft interaction 
CTsoft is considered nonperturbative and thus noncalcula- 
ble, the jet production cross section djot is assumed to be 



given by perturbative QCD (pQCD) for transverse mo- 
mentum transfer px > Po- The two parameters, CTsoft and 
Po, are determined phenomenologically by fitting the ex- 
perimental data of total p+p(p) cross sections within the 
two-component model [|lO|- [l3| , |l6| , p7t . 

The cut-off scale po, separating nonperturbative and 
pQCD components, could in principle depend on both en- 
ergy and nuclear size. Using Duke-Owens parameteriza- 
tion psf of parton distributions in nucleons, it was found 
in the HIJING jl^ model that an energy independent 
cut-off scale po = 2 GeV/c is sufficient to reproduce the 
experimental data on total cross sections and the hadron 
multiplicity in p + p{p) collisions, assuming that the soft 
cross section (Xsoft is also constant. Since then, analysis 
of recent experimental data from deep-inelastic scatter- 
ing (DIS) of lepton and nuclcon at HERA indicated |l^] 
that gluon distribution inside a nucleon is much larger 
than the DO parameterization at small x. Many new pa- 
rameterizations of the parton distributions have become 
available. Using the Gluck-Reya-Vogt (GRV) parameter- 
ization |ig| ] of parton distributions and following the same 
procedure as in the original HIJING jl^, we find that 
one can no longer fit the experimental p + p{p) data us- 
ing a constant cut-off scale po within the two-component 
model. One has to assume an energy dependent cut-off 
scale po(-\/s)- Because of the rapid increase of gluon dis- 
tribution at small x, we find that the cut-off po(\/s) ^las 
to increase slightly with energy in order to fit the exper- 
imental data. 

Shown in Fig. ^ is the calculated central rapidity den- 
sity. 
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for p + p(j)) collisions as a function of energy y^s, where 
(n)s — 1.6 and (n)h — 2.2 represent particle production 
from soft interaction and jet hadronization, respectively. 
The jet cross section in lowest order of pQCD is given by 
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CTjet = K dprpdyidy2 
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a,b,c,d 



XlX2fa(.Xi)Jb{x2) 



dcFa. 



h — ^cd 



dt 



(2) 



with the GRV parameterization of parton distribu- 
tions fa{x), a factor of 2 and an energy-dependent 
cut-off scale 

Po(Vs) = 3.91 - 3.341og(log Vs) + 0.98 log^ (log \/s) 

0.23 log^ (log Vi). (3) 

Assuming eikonalization of hard and soft processes, the 
total inelastic j3+p(p) cross section in this two-component 
model is 1T^, 



= d'b 



2 _ g-(o'jot-Ho'Boft)T'jviv(f)) 



(4) 



where T]\[i^{b) is the nuclcon-nucleon overlap function 
and tJsoft =57 mb represents the inclusive soft cross sec- 
tion. Notice that the energy-dependence of po{\/s) is 
quite weak, ranging from po = 1.7 GeV/c at ^/s — 20 
GeV to 3.5 GeV/c at Vs = 5 TeV. 
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FIG. 1. Charged particle rapidity density per participat- 
ing nucleon pair versus the cm. energy. The RHIC data 
(filled circle and up-triangle) for the 6% most central Au+Au 
are compared to pp and pp data (open symbols) |2c|-^ and 
the NA49 Pfe -h P6(central 5%) data H (filled square). The 
two-component mini-jet model with and without shadowing 
is also shown. The shaded area for central Au + Au colli- 
sions corresponds to the range of gluon shadowing parameter 
Sg = 0.24-0.28 [Eq. (||)]. 

To extrapolate the two-component model to nuclear 
collisions, one assumes that multiple mini-jet production 
is incoherent and thus is proportional to the number of 
binary collisions iVbinary The soft interaction is however 
coherent and proportional to the number of participant 
nucleons A'pait according to the wounded nucleon model 
. Assuming no final state effects on multiplicity from 
jet hadronization, the rapidity density of hadron multi- 
plicity in heavy-ion collisions is then, 



= ^(^part)('^)s + (n)/i(A^bmary) 



(5) 



where CTj^^(s) is the averaged inclusive jet cross sec- 
tion per NN in AA collisions. The average number of 
participant nucleons and number of binary collisions for 
given impact-parameters can be estimated using HIJING 
Monte Carlo simulation. If one assumes that the jet pro- 
duction cross section a^^{s) is the same a.smp + p colli- 
sions, the resultant energy dependence of the multiplicity 
density in central nuclear collisions is much stronger than 
the RHIC data as shown in Fig. ||. Therefore, one has 
to consider nuclear effects of jet production in heavy-ion 
collisions. 

In high-energy nuclear collisions, multiple mini-jet pro- 
duction can occur within the same transverse area. If 
there are more than one pair of mini-jet production 
within the transverse area given by the jet's intrinsic size 
tt/pq, jet production within this area might not be in- 
dependent any more [l^ . If such a criteria is used for 
independent jet production within one unit of rapidity, 
one can then obtain a cut-off scale po in a so-called final 
state saturation model |25|, 
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(6) 



that also depends on nuclear size for central heavy-ion 
collisions. This cut-off scale, though increasing with nu- 
clear size, ranges from 0.7 GeV/c at ^/s — 20 GeV to 
2.2 GeV/c at ^/s = 5 TeV for central Au + Au collisions, 
which is much smaller than what we have obtained in 
Eq. ^ by fitting p ~\- p{p) data. Therefore, if we apply 
the two-component model to heavy-ion collisions with the 
same cut-off scale in Eq. (||) as determined in p + p{p) 
collisions, the criteria for independent jet production will 
never be violated. Instead, we will assume the cut-off 
scale to be independent of nuclear size in this paper. 

In principle, jets produced in the early stage of heavy- 
ion collisions will also suffer final state interaction and 
induced gluon bremsstrahlung. For an energetic jet, this 
will lead to induced parton energy loss ||2^,|2^ and the 
suppression of large transverse momentum hadrons [ p8[ . 
Such a jet quenching effect could also lead to increased 
total hadron multiplicity due to the soft gluons from 
the bremsstrahlung. However, a recent study [p9| of par- 
ton energy loss in a thermal environment found that the 
effective energy loss is significantly reduced for less ener- 
getic partons due to detailed balance by thermal absorp- 
tion. Thus, only large energy jets lose significant energy 
via gluon bremsstrahlung. Since the production rates of 
these large energy jets are very small at the RHIC en- 
ergy, their contributions to the total hadron multiplicity 
via jet quenching should also be small. Similarly we also 
assume that parton thermalization during the early stage 
contributes little to the final hadron multiplicity. 

One important nuclear effect we have to consider in 
our two-component model is the nuclear shadowing of 
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parton distributions or the depletion of effective parton 
distributions in nuclei at small x. Such a nuclear shad- 
owing effect in jet production can be taken into account 
by assuming modified parton distributions in nuclei, 



(7) 



Using the experimental data from DIS off nuclear tar- 
gets and unmodified DGLAP evolution equations, one 
can parameterize R^{x^ Q^) for different partons and nu- 
clei 1^,^. Recent new data however indicate that 
the simple parameterization for nuclear shadowing used 
in HIJING jl^ is too strong for heavy nuclei. In this 
paper, we will use the following new parameterization. 



R^{x) = 1.0 + 1. 19 log^/*k(a;3 



1.2x^ + 0.21a;) 



- Sg (^1/3 _ 1)0.6(1 _ 3 5^) exp(-2;V0.01) (8) 

with Sq = 0.1 for all quark distributions as shown in 
Fig. H (solid lines) in comparison with the experimen- 
tal data |30|. Also shown in the figure are parameter- 
izations (dashed lines) by Hirai, Kumano and Miyama 
(HKM) and the old HIJING parameterization 
The shadowing in the old HIJING parameterization (dot- 
dashed lines) is apparently too strong for heavy nuclei. 
The HKM parameterizations also take into account con- 
straints by momentum sum rules, as similarly in the orig- 
inal parameterizations by Eskola,Kolhinen and Salgado 
(EKS) [0. For the purpose of this paper, one can ne- 
glect the scale dependence of the shadowing. 
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FIG. 2. Ratio of nuclear structure functions as measured 
in DIS. Solid lines are the new HIJING parameterization 
[Eq. 1^ , dashed lines are the HKM parameterization ^] and 
dot-dashed hnes are the old HIJING parameterization |^^. 
The data are from Ref. po| . 

The nuclear shadowing for gluons is somewhat con- 
strained by the momentum sum rules in the HKM param- 
eterization. However, the constraint is not very strong. 



leaving a lot of room for large variation of gluon shad- 
owing. Shown in Fig. ^ are the shadowing factors for 
gluon distribution from EKS and HKM parameteriza- 
tions. They both have strong anti-shadowing around 
X ^ 0.1. The stronger anti-shadowing in EKS param- 
eterization is due to additional constraints by the 
dependence of F2{Sn)/ F2{C), assuming the same un- 
modified DGLAP evolution equation for parton distri- 
butions of a nucleon. Since gluon-gluon scattering domi- 
nate the jet production cross section, such a strong gluon 
anti-shadowing leads to larger jet cross section and thus 
larger hadron multiplicity than in the case of no shad- 
owing at the RHIC energies. Such a scenario within the 
two-component model is clearly inconsistent with the ex- 
perimental data. We therefore propose the following pa- 
rameterization for gluon shadowing. 



Rf{x) = 1.0 + 1.191og^/*k {x^ 



(A 



1/3 



1)0-6(1 



1.2x^ -H0.21j;) 
1.5a;°-3^)exp(-a;V0-004), (9) 



with Sg = 0.24-0.28. This is shown in Fig. || as the 
solid lines. The hadron multiplicity density in the two- 
component model using the above gluon shadowing is 
shown in Fig. |l|. The shaded area corresponds to the 
variation of Sg = 0.24-0.28. The RHIC data thus in- 
dicate that such a strong gluon shadowing is required 
within the two-component model. If one assumes the 
same gluon shadowing as the quarks in Eq. (^), the re- 
sultant dN/drj is only slightly smaller than the one with- 
out shadowing. Such a contraint on gluon shadowing is 
indirect and model dependent. It is important to study 
directly the gluon shadowing in other processes in AA or 
pA collisions. 




10 X 10 I 

FIG. 3. Ratios of gluon distributions in different nuclei 
given by the new HIJING [Eq. ^] (solid line, the shaded 
area corresponds to Sg = 0.24-0.28), old HIJING (l|] 
(dot-dashed), HKM Q (dashed) and EKS [|l) (dotted) pa- 
rameterization. 
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To take into account the impact-parameter dependence 
of the shadowing, we simply replace the shadowing pa- 
rameters Sa in Eqs. (||) and by 



Sa{b) 



(10) 



where Ra = 1.12A^^^ is the nuclear size. With this 
impact-parameter dependence, the calculated jet cross 
section cr^^(s)/CTi„ will also depend on the centrality of 
heavy- ion collisions, decreasing from peripheral to central 
collisions. One can then calculate the centrality depen- 
dence of the hadron multiplicity density. The results are 
shown in Fig. ^ as a function of A^part for Au + Au col- 
lisions at y/s = 56, 130 and 200 GeV. The shaded areas 
again correspond to the variation of the gluon shadowing 
parameter Sg = 0.24-0.28. Within statistical and sys- 
tematic errors, the two-component mini-jet model with 
impact-parameter dependent parton shadowing describes 
the PHOBOS and PHENIX data |,| at ^/s = 130 GcV 
well. 
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FIG. 4. The charged hadron central rapidity density per 
participant nucleon pair as a function of the averaged num- 
ber of participants from the two-component model (shaded 
lines) , two-parameter fit [Eq. (|ll])] (dot-dashed lines) and par- 
ton saturation model bl as compared to experimental data 

To illustrate the effect of the impact-parameter depen- 
dence of the parton shadowing, we compare the results 
with a two-component parameterization. 



dNch 
dr] 



1 



(11) 



shown as dot-dashed lines, where the two parameters, rig 
and n/i, fixed at each energy by values of dN^/dij in 
p + p collisions and the most central Au + Au collisions, 
are assumed to be independent of the centrality. The 
increase of 2dNch/ drj / {Np^^rt) with {Nps,H) is driven only 
by the centrality dependence of (A'^binary) / (iVpart) in this 
two- parameter fit. Comparing to such a two-parameter 



fit, the two-component mini-jet model has a flatter cen- 
trality dependence at high energies because the effective 
jet cross section decreases from peripheral to central colli- 
sions due to the impact-parameter dependence of parton 
shadowing. The better agreement between the experi- 
mental data and the two-component mini-jet model at 
= 130 GeV is another indication of strong nuclear 
shadowing of the gluon distribution in mini-jet produc- 
tion. 

Similar centrality dependencies are also predicted by 
other models p5| , p3| , in particular the initial-state par- 
ton saturation model 1^,1). It is based on the nonlin- 
ear Yang-Mills field dynamics |l^j3^ assuming that non- 
linear gluon interaction below a saturation scale ~ 
asxGA{x,Q^)/TTR\ leads to a classical behavior of the 
gluonic field inside a large nucleus, where Ga{x,QI) is 



the gluon distribution at a; = 2Qs/\fs. Assuming parti- 
cle production in high-energy heavy-ion collisions is dom- 
inated by gluon production from the classical gluon field, 
one has a simple form ||^ for the charged hadron rapidity 
density at ?7 = 0, 



dNr 



part/ 



di] 



So 



A, 
+ 2 log 



A/2 
S 

So 



log 



Qos 
^QCD 



(12) 



• PHOBOS Au-l-Au 



with c w 0.82 This is shown in Fig. ^ as solid lines. 
Here, Aqcd = 0.2 GeV, A = 0.25 and the centrality de- 
pendence of the saturation scale Qq^ at y/so = 130 GeV 
is taken from Ref. M. 
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FIG. 5. The ratios of charged hadron multiplicity density 
in Au + Au collisions at different energies as a function of the 
averaged number of participants shown with PHOBOS data 



Comparing the two model results in Fig. U, one no- 
tices that the saturation and two-component model agree 
with each other in most regions of centrality except very 
peripheral and very central collisions. In central col- 
lisions, results of saturation model tend to be flatter 
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than the two-component model. In this region, there 
are still strong fluctuations in parton production in the 
two-component model through the fluctuation of A^binary 
while A'part is limited by its maximum value of 2A. 



That is why dNch/dr]/{N^ 



part , 



continues to increase with 



(A^part) in the central region. Such a fluctuation is not 
currently taken into account in the saturation model cal- 
culation. More accurate measurements with small errors 
(less than 5%) will help to distinguish these two different 
behaviors. For peripheral collisions, saturation model re- 
sults fall off more rapidly than the mini-jet results. How- 
ever, the experimental errors are very big in this region 
because of large uncertainties related to the determina- 
tion of the number of participants. Therefore, it will be 
very useful to have light-ion collisions at the same energy 
to map out the nuclear dependence of the hadron multi- 
plicity in this region. An alternative is to study the ratios 
of hadron multiplicity of heavy-ion collisions at two dif- 
ferent energies as a function of centrality. In this case, the 
errors associated with the determination of centrality will 
mostly cancel. Shown in Fig. ^ are the ratios of hadron 
multiplicity at three different energies as a function of 
the averaged number of participants predicted by satu- 
ration and two-component model. One notices that while 
the results from saturation model have the same central- 
ity dependence at all three energies the two-component 
model predicts slightly different behavior at different en- 
ergies, indicating the energy dependence of the mini-jet 
component. So the ratios given by the saturation model 
are almost independent of centrality. On the other hand, 
two-component model predicts noticeable centrality de- 
pendence of the ratios. This is especially true for the 
ratio between colhsions at ^/s = 200 and 56 GeV. 

It is interesting to point out that in the saturation 
model that assumes a particle production mechanism 
dominated by coherent mini-jet prodctuon below the sat- 
uration scale Qs , the value of Qs determined in Ref. [HJ|] 
is much smaller than the cut-off po in the two-component 
model constrained by the p+p{p) data. As demonstrated 
in this paper, the number of mini-jet production below 
such scale is still very large and should contribute to the 
final hadron multiplicity. 

In summary, we have studied the energy and centrality 
dependence of the central rapidity density of hadron mul- 
tiplicity in heavy-ion collisions at RHIC energies within 
a two-component mini-jet model. As a consequence of 
the latest parameterization of parton distributions 
which have a higher gluon density than the old param- 
eterization |l8) used in previous studies the cut-off 
scale that separates soft and hard processes is found to 
increase slightly with energy in order to fit the p + p{p) 
data. The cut-off scale, however, is still large enough that 
the independent jet production picture is still valid. With 
a new parameterization of nuclear shadowing of parton 
distributions in nuclei, we also found that RHIC data 
require a strong shadowing of gluon distribution. Using 
this strong gluon shadowing with an assumed impact- 
parameter dependence, the predicted centrality depen- 



dence of the hadron multiplicity agrees well with the re- 
cent RHIC results. We have also compared our results 
with the parton saturation model |^,^ . We point out that 
in order to differentiate the two models one needs more 
accurate experimental data in both the most central and 
peripheral regions of centrality or study the centrality 
dependence of the ratios at different colliding energies. 
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